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Abstract 

This is the third paper in a series outlining an algorithm to con- 
struct finite states of quantum gravity in Ashtekar variables. In this 
paper we treat the case of the Klein-Gordon field quantized with grav- 
ity on the same footing. We address the full theory, outlining the 
solution to the constraints and the construction of the corresponding 
wavefunction of the universe. The basic method for the full theory is 
to expand the constraints relative to the solution for the pure Kodama 
state and rewrite them in the form of a generalized nonlinear group 
transformation of the CDJ matrix, viewed as a nine-dimensional vec- 
tor. We then outline a prescription for finding the fixed point of the 
flow, and the corresponding generalized Kodama state for the full the- 
ory is constructed. The flnal solution is expressed in an asymptotic 
series in powers of model-specific matter inputs, suppressed by a small 
dimensionless constant, relative to the pure Kodama state. We discuss 
this expansion from different perspectives. Lastly, we explicitly show 
how the solution to the quantized constraints establishes a wavefunc- 
tion of the universe with a predetermined semiclassical limit built in 
as a boundary condition on quantized gravity in the full theory. 
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1 Introduction 



A linearized treatment of general relativity in metric variables involves an 
expansion of the spacetime metric in fluctuations about a specified back- 
ground metric. For a flat Minkowski background, a ground state of metric 
relativity, the fluctuations have the interpretation of gravitons h^i, , given by 

9ixu =ViJ.u + y^hij,^- (1) 

When (1) is substituted into the action for general relativity, an infinite 
series is obtained in powers of the coupling constant VG. It was first pointed 
out by Pauli that a loop expansion of the effective action of metric general 
relativity would be nonrenormalizable due to the negative mass dimension 
of VG. In this paper we will try to make sense of an analogous expansion 
by treating the CD J matrix as a metric in S'C/(2)_ via the expansion 

6 -A. 

^ae = -j{Sae + -^eae). (2) 

In (2) 

Gae is the CDJ deviation matrix, which parametrizes departure from 
the CDJ matrix for the pure Kodama state 6ae- One difference between (2) 
and (1) is the mass dimension of the coupling constants ([A] =2, [G] = —2). 
We will ultimately suggest the hypothesis for a renormalizable description 
in Ashtekar variables due to the positive mass dimension of A. 

It was argued in [1] that a linearization of the pure Kodama state about 
DeSitter spacetime is unnormalizable in the inner product of the theory of 
linearized gravitons in the Lorentzian theory. However, it was also shown 
that the pure Kodama state ^xod is the semiclassical manifestation of De- 
Sitter spacetime [2], [3]. This begs the question, when one wishes to assess 
the suitability of the pure Kodama state ^Kod as a vacuum state of general 
relativity with a good semiclassical limit, as to whether one should rather 
expand fiuctuations relative to '^Kod rather than expand ^Kod itself. 

In [4] Witten argues that "^Kod cannot serve as the ground state for quan- 
tum gravity by analogy with Yang-Mills theory, the Chern-Simons state of 
which contains severe pathologies such as unnormalizability and nonunitar- 
ity. However, Randono presents arguments in [5], [6] which addresses these 
issues when one interprets '^Kod essentially in analogy to plane wave momen- 
tum eigenstates. Furthermore, Smolin and Soo in [7] present the hypothesis 
for ^Kod as the ground state of qiiantum gravity with A term by studying 
the behavior of small perturbations of ^Kod due to coupling of the theory 
to matter. In this work the perturbations were measured relative to the 
self-dual sector for DeSitter spacetime, parametrized by 
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A 

It was shown for small perturbations that the matter degrees of freedom 
roughly correlate to g"* lending some credence to the notion of ^Kod as 
a proper ground state. When one applies the (nondegenerate) Ashtekar 
magnetic field to (2) and uses the CD J Ansatz [11], [13], one obtains 



al + jBl = -eaeBl. (4) 

Equation (4) shows that the CDJ deviation matrix eae is simply g"' with 

the spatial index i projected into SU{2). As pointed out in [7] it would be 
necessary to study the higher order corrections to the semiclassical limit in 
order to further investigate this hypothesis, but no method exists to do so. 

The ultimate purpose of this paper is to illustrate a new procedure by 
which the constraints of general relativity can be solved in the full theory, 
which enables the construction of the generalized Kodama state '^oKod [H] 
devoid of ultraviolet singularities. The claim is that ^cKod amounts to an 
expansion about '^Kod in terms of the matter content of the model. Our 
new method will illustrate a systematic algorithm to obtain this expansion 
to any order desired, which should enable the hypothesis in presented in [7] 
to be tested for a large class of model-specific matter couplings. 

We will first rewrite the nine equations stemming from the quantum 
constraints in terms of the deviation matrix e^fc) using the results from [11]. 
Section 2 puts the constraints into standard form amenable to the expansion 
(2), which is carried out in sections 3 and 4. We then create a library of 
terms which can be utilized to examine more general models of interest. 

In section 5 we represent the constraints as a linear operator acting on 
eae plus nonlinear corrections suppressed by a small dimensionless coupling 
constant GA. In a matrix representation acting on a nine dimensional space 
of deviation matrix elements eae, this has the interpretation of the flow of 
model-specific matter sources under the action of a new type of nonlinear 
transformation. The linear part of the transformation has the representation 
of a nine by nine matrix of functional integro-differential operators. The 
differentiation occurs in the direction of spatial position x e T, and also in 
the direction T of the functional space of quantum fields Af . It is convenient 
to think of these variables as comprising the base space of a fibre bundle, 
with eae as a section on the bundle undergoing a nonlinear transformation. 

The solution to the quantum constraints then amounts to finding the 
fixed points of this transformation, using the matter contributions as a driv- 
ing input. A full solution requires inversion of the linear part of the transfor- 
mation followed by iterative perturbation expansion. In section 6 we suggest 
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a method to carry out this inversion, in analogy to the solution of the Gauss' 
law constraint discussed in [11]. 

The precise determination of the propagation kernels required for full 
inversion require solution of linear first-order and a linear second-order dif- 
ferential equations with nonconstant coefficients in the full theory, which 
we do not carry out in this paper. In section 7 we illustrate the solution 
for the deviation matrix by nonlinear iteration of the aforementioned linear 
transformation. We also illustrate the construction of the generalized Ko- 
dama state for the full theory for gravity quantized with the Klein-Gordon 
scalar field, as an asymptotic expansion relative to the pure Kodama state 
in powers of GA. The coefficients of the expansion involve propagation of 
the nine-dimensional 'generalized' matter charges, both in position space S 
and also in the direction of the functional space of fields Fa- 



2 Standard form for qq, qi and q2 in Ashtekar con- 
nection variables 

As a review from [11], the coefficients of singularity determined from the 
quantum Hamiltonian constraint are given by 



qo = detBlyar"^ + Adet*) + 0^0 = 0; 

k d d 



Q2 



6 dAf dA^, ' 

' 3 



ei^k^^'^'^ceBl) + 36 = 0. (5) 



It will be useful to express the constraints (5) in the form of an operator 
acting on the nine-dimensional vector space of CD J matrix elements. This 
requires that the constraints be put into standard form. The kinematic 
constraints, being linear in the CDJ matrix, are already in standard form. 
This leaves the terms of the quantum Hamiltonian constraint qq, q\ and q2. 



2.1 Semiclassical term go 

Starting with the semiclassical term Hd = qo, which is of zeroth order in 
singularity, we have 
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qo = detB {Var-^ + Adet*) + G^o • (6) 

Note that the condition qo = imphes a well-defined solution only for nonde- 
generate magnetic fields when there is matter in the theory (e.g. Clo ^0). 
The case of degenerate B is not applicable to the generalized Kodama states. 



2.2 First-order singuleirity qi 

The first-order singularity term, given by qi is a first-order nonlinear differ- 
ential equation given by 

+^:^i'^ijke''"'^ce^bfB',Bj) + GQi. (7) 

Equation (7) can be thought of as the functional divergence of linear and 
quadratic CD J matrix terms with respect to 2-index and 4-indcx SU{2)- 
tensorial indices, which resembles a minisuperspace equation but is actually 
for the full theory. The matter contribution J7i is present for matter models 
with Hamiltonians quadratic in conjugate momenta. The functional form 
of the solution is the same for all x (dependence of which is suppressed), 
since the x dependence must be frozen in order that the integration can be 
carried out in the direction of the functional space of fields A", ^. We now 
place the gravitational part into standard form. 

2.3 Linear parts of qi 

Expanding first the part of (7) linear in CDJ matrix elements, we have 

{qi)linear = ^(eij^e'^'^BiB^.^.e) + eij^e'^'^D^Ji^aeBl) 
= e^j^e^^^^iDiiB', + BlD^J,)^,, + e,,fc6"^^Z?S5^*„e + e.^^e^^'^BiBl^^,, (8) 
Starting with the first term on the right hand side of (8), 

= (e,,fc6^-^')(6'^^'=6e„d)Afi3c'^5e = (- 2^) (- 1) (^^^^^ " S^St)AlB^^,, 

= 2{6l6^, - 6'^6:;)C^^,e = 2(<5^rC - C^)^,^- (9) 
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Moving on to the second term on the right hand side of (8), 



= {eijke''%^''e^d)BiAf = {26\){25\)BiAf = ACl^^ (10) 
Moving on to the third term on the right hand side of (8), 

= {eijke''''W''ecM)AfBi^ae = {-2Si){-26^a)AfBi^ae = ^C^^ae = 4C,^^'6e(ll) 
Finally, the fourth term on the right hand side of (8) 



e.j.e'^'^BiB'^-^^,, = idetB)e'^%UB-')f^^be (12) 



So the total contribution to qi linear in the CDJ matrix elements is given 
by the sum of (9), (10), (11) and (12), or 

{qihnear = ^ (€.,fc€'^'"=S^i?,^*fee) + 6.,fce'^'"=Dg (*aei?:) 

= {detB)e'^'%^{B-')f^ + (2(<5^rC - Cl)^,,) + (4C,^*6e) + i^C'e^be) 



{detB)e"''^e,^{B-^)f-^ + 6C', + 26',tiC ^,,.{^3) 



2.4 Quadratic parts of gi 

Moving on to the quadratic CDJ matrix contributions to qi, we expand the 
part of (7) quadratic in CDJ matrix elements. 



Ad ■ 

{qi)quad = ^^]4^(^«ife^" "-^ce'^bfB^Bjj 



I 



= ^e,j,e<^'^{D':iBj + B^,D%)^,,^,f + ^^,.,^"''^5,^5}^*,^*./ (14) 
Starting with the first term on the right hand side of (14), 

= ^{eijke''%''''eead)AfBj^S^^S,f = ^i2S^^){-l){S% - S',S^)AfBj^,,^,f 

= ^i-2){5l5^, - 5'^6^)Cp,,^,f = j(-2)(5,^q - 5,^Cj)*ee*v(15) 
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Moving on the the second term on the right hand side of (14), 



id rik-i 



And finally, the last term on the right hand side of (14), 



So the total contribution to qi quadratic in the CDJ matrix elements is given 
by the sum of (15), (16) and (17), or 



A d 



A 



d 



{Qi),uad = --^ieijke''''^ce^,fB',Bj) = -(detB)e^'-efed{B-')t-^^ce^,f 

+ (^^i-2){6lCf - <5,^C))*ee*v) + (^(2)(4Ce " S}Ct)^ce^hf 



(18) 



It can be shown by reshuffling indices that the last two terms on the right 
hand side of (18) are equal, leading to the result that 



A 



{qi),uad= ^(detB)6«''V<i(^"')'^+A(5^C,^-<^^Ce')J*ce*6/ (19) 



2.5 Total contribution to functional divergence terms 

So the total contribution to the term first-order in singularity is given by 
the sum of (13) and (19) 



Qi 



(det5)e«^^edec(S 



-l\d 



d 
dA§ 



+ + 25lixC 



he 



+^[AeiB)e'^''^eafe{B-^)i^ + K^)Cl - 5^C,^)] ^ee^v + GQi (20) 
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2.6 Functional Laplacian term 

The functional Laplacian term q2 is given by 



92 



A_d d_ 

6dAfdA'i 



(€,jfce"^^*eei?e ) + 36 = A^^^e + 36 



be. 



(21) 



This can be expanded as 



q2 - ei^ke^^^Bl-^^^!^ + eijke'""'D';l^^ce 



d d 



abc T^ki 



d 



d d 



dAf dA]' 



(22) 



(I) We now evaluate the contributions to q2 first-order in derivatives. Start- 
ing with the second term of the right hand side of (21), 



-2K«-W(^^); = -2[( 



which is traceless. The third terms on the right hand side of (21) and the 
second term are both equal. 

Expanding the third term on the right hand side of (21), 



d 



d 



which is equal to (23). 

(II) We now look at the contribution to q2 of second order in derivatives. 
The last term of (21) is given by 



d \c 



d 



' 3 



abc 



d_ d_ 

dA^dM 



{e'^'^'diAem + ^e'='-eed/^f^4),(25) 



where we have made use of the unconventional calculus, first used in [11], 
that functional and spatial differentiation must commute. Continuing on, 
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abc 



d_ d_ 



{e^^--diA,m + \e^''^e,afAfAi:) 



= le,,,e'^'"^e''^ee,f^{5t5jAi + Afsld^) 



abc klm 



edfisisisidi+disididi,) 



= ie,,.,€'^'"=e*=^\e6a + le,,ke''''e'^''e,ab = (-3)(-25,^) + 3(25,^) = 125,^ (26) 
which gives a total, upon combining (23), (24), (25), (26), of 



A 



92 



*ce + 36, (27) 



3 Expansion of the quantized kinematic constraints 
relative to the pure Kodama state 

Now that the coefficients of the singularities for the quantum Hamiltonian 

constraint have been put into standard form, we next expand these coeffi- 
cients relative to the pure Kodama state ^Kod- The parametrization of the 
CDJ matrix consists of the expansion 

*ae = -(«~^<^ae+eae) (28) 

where in (28) we have defined k = (1/6) A^, and where the deviation matrix 
€ae consists of nine independent elements. The elements eae arc independent 
of each other until a solution to the quantum constraints is found, whereupon 
its degrees of freedom become saturated. 

Let us begin with the kinematic constraints, which are simultaneously 
satisfied at the semiclassical and at the quantum levels (SQC) due to be- 
ing linear in conjugate momenta [12], and hence linear in the CDJ matrix. 
Note that the mass dimensions are [k ] — [^ae] — — 

2. Starting with the 

quantized diffeomorphism constraint [11], 



^dae^ae = —GTQd- (29) 

'^It will be convenient to use k instead of A in many of the computations that follow, 
to avoid carrying around various numerical factors of 6. 
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The antisymmetric elements of the CD J deviation matrix Eae are determined 
by the local matter momentum Hi. Since the constraint must be identically 
satisfied at all points x in the spatial hypersurface S^, labeled by time t, 
these elements are uniquely globally determined in the functional depen- 
dence upon the fields {Af{x), <^"(.x)) indicated at all points. Note that since 
this is a linear constraint it can be solved explicitly when the Ashtekar mag- 
netic field is nondegenerate (invertible) . This invertibility condition is a 
direct consequence of the presence of local matter fields in the theory [13]. 
Also, the pure Kodama part of the equation has vanished. A final note for 
checks of consistency in the last line of (29) is that the mass dimension of 
Hi is [Hi] = [0] -I- [5] -I- [tt] = 1 -I- 1 -I- 2 = 4, causing the mass dimensions to 
balance. 

Now for the quantized Gauss' law constraint, 

BiDi^ae + GQa = B^Di [-K'^Sae - e„e] + GQ^ = 
BiDiSae + BiDiCae = GQ„ ^ + fi^ACae = GQa (30) 

Note again that the pure Kodama contribution has dropped out, leaving 
the remaining parts e^e to be determined by the matter fields. This leads 
to three first-order differential equations in nine unknowns [11] 

BiDieae = [Safg^ + C/^) 6/, = = GQa (31) 

where C = = A^B^ is a SU (2) ® SU (2)-valued object which acts as a 
connection for parallel transport along the directions t^, where e = 1,2, 3 
represent SU(2) indices. The integration of this partial differential equation 
requires as boundary conditions the specification of the values of three ele- 
ments Cae on three linearly-independent 2-dimensional hypersurfaces dual to 
the direction of path-ordered integration of the partial differential equation 
(31) and the value of all elements not being solved for via this constraint 
(six) everywhere in S. Three elements of eae are already fixed by (29). So 
let us decompose the CDJ matrix into its symmetric and its antisymmetric 
parts 

^ae = V'ae + CaedV'd (32) 

Then Gauss' law constraint can be solved as outlined in [11] for three of the 
symmetric elements ^ae in terms of the antisymmetric elements ip^^ which 
are uniquely determined by the local matter momentum, and the remaining 
three symmetric elements. This yields 

BiDi{'^ae + WV'd) = GQa 

BlD,^ae = GQa - GB^Di (w^^) • (33) 
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So the Gauss' law constraint has been reduced via (33) from nine to three 
equations in six unknowns. These equations can be solved explicitly for 
three elements in terms of the remaining three, which are determined by the 
Hamiltonian constraint. Lastly, regarding mass dimensions, the parameters 
of path-ordered integration have [t^ = [d] + [B] = 1 + 2 = 3 and the SU{2)_ 
charge has [Q] = [(/)] + [ir] = 1 + 2 = 3, causing the mass dimensions on 
both sides of (30) to balance. In this work we will not use the diagonal CDJ 
matrix elements to solve the Gauss' law constraint as shown in [11], but 
rather the upper off-diagonal elements *i2, *23, ^3i- 



4 Expansion of the Quantized Hamiltonian con- 
straint relative to the pure Kodama state 

We will now expand the Hamiltonian constraint in terms of the Ansatz 
(2). This is not a linearization, since we will carry out the expansion to 
all orders. The polynomial nature of the constraints in Ashtekar variables 
ensures termination of the series expansion in A unlike in metric variables 
where the corresponding series is infinite. The standard form of the Hamil- 
tonian contribution to the constraints is then given, rewriting (48) and (18) 
for completeness, by 

qo = detB (Far* + Adet^') + GJ^o (34) 



Qi 



(detB)e«^^edec(5-^)f ^ + 6C', + 2<5^trC 



be 



A 

q2 = -p: 



*ee + 36, (36) 



4.1 Expansion of q2 relative to the pure Kodama state 

We now expand the CDJ matrix about the isotropic pure Kodama solution 
in the quantized Hamiltonian constraint. 
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A 

(12 = 

6 



^abc r) 



'dA)dAl 



dAl 



+ I25a 



(--^(^ae - eae) + 36 



= ^(-|<5„e)(12<5ae) + 36 + ^A^eEae = ^A„eeae(37) 

The pure Kodama contribution has cancelled out the inhomogeneous factor 
of 36. Equation (37) can be thought of as a set of nine individual 'gen- 
eralized' functional Laplacian operators A^e, each acting on an individual 
element of e^e- The functional Laplacians are given by 



e^^ + mAl-S^At)^ 



dA)dAl 



dAl 



(38) 



Since each elements of A^e is at worst a linear record-order differential op- 
erator at each point x with respect to the functional space of fields, its 
inversion can be carried out in principle by heat kernel methods and ex- 
istence theorems regarding the solutions can be applied. We relegate this 
technical step of the full superspace theory to future publication. 



4.2 Expansion of the nonderivative parts of qi relative to the 
pure Kodama state 

We are now ready to carry out the expansion on the nonderivative parts of 

^ae = -{-5ae + eae) (39) 
K 

where n = (1/6) A. Starting with the nonderivative linear term in giin the 
CDJ matrix, 

(6C) + 25)ixC)^ae = (6C) + 26)iiC){-\f - e^f) 

= --{GC^fSbf + 2S'}Si,ftTC) - (6C) + 2S'}tTC)€bf 

= -(^)trC- {pC) + 25)tTC)€bf. (40) 
Now the for the nonderivative part of qi quadratic in the CDJ matrix, 

= QKCff{-5ae + eae) {-Sbf + ^bf) 
= GnC^f (^SaeSbf + -Saeebf + -^aeSbf + eaeCft/) • (41) 
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Concentrating on the first term of (41), the contribution due to the pure 
Kodama state is given by 



= ^(3trC-<^.,C,^) = l^ (42) 

Note that the factor of 12 on the last Hne of (42) is precisely that needed to 
cancel the leading order term of (40). 

Next we move on to the terms in (41) linear in 6(16 ■ Contrary to the 
inclination to double this term, thus avoiding its explicit calculation, let us 
proceed with its computation. 



= 6(3C} - C'})ebf = 12tr(Ce) (43) 
Moving on to the next linear term of (41), 



hi 

= 6((trC)tre - eabS^Cpbf) = 6(trC)tre - 6tr(Ce) (44) 

(44) is not quite the same as the previous term (43), which may be counter 
to any preconceived intuition on symmetries. This is a manifestation of a 
property which has a statistical interpretation, known as skewness. 
Lastly, there is a term of (41) quadratic in e 

G^Cffeaeebf = {S^C'f - 5lC1)e^ebf = AC^e^eeft/. (45) 

This is the first term featuring a departure from linearity. Notice that 
this term is suppressed by a factor of order ~ A which is small, therefore 
one would expect a linearized treatment of this particular term (45) to be 
extremely accurate. 

To recapitulate, combining all contributions, the nonderivative part of 
the expansion of qi is given, combining (40), (42), (44) and (45), by 

+ + 12tr(Ce) + 6(tre)(trC) - 6tr(eC) + 6AcCj^e„ee6/ 



= -(eq + 25)tiC)ebf + 6(tre)(trC) + 6tr(eC) + &KC%taee 



bf 



6tr(eC) + 6(tre)(trC) - {QC) + 25}trC)eb/ + KC%eaeebf, (46) 
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where we have used the cychc property of the trace. 



4.3 Expansion of the functional derivative parts of qi relative 
to the pure Kodama state 

Under the Ansatz 

^bf = -{\f+e,f) (47) 

rv 

and using (48), we must first focus on the gravitational part of 



A 



d 



+ :i{detB)e^'-efUB-')f-^ + A{6}C^^ - S}C',)^ce^,f\ + GJ^i.(48) 



This gives 



e«^'=(detS)(S-l)f ^ [eecd^be + ^l^efed^ce^bf 



(49) 



Starting with the hnear term of (49), we have 



,a6c 



{detB){B 



-l\d 



d 



dA' 



a^ecdi ebe) = -(detB)e"''^eecd(5 



-l\d 



{detB)e^'%dciB-')f^e,e = (det5)(5,»5^ - S^,5l){B~')f 



d 

d 



(detB) {B 



-l\b 



dA, 
d 



(■he 



e^. (50) 



Moving on to the quadratic term of (49), we have 



^/.6«^^(det5)(5-i)f^6/ed*ce*;./ = ^«(detS)e"^^6^ed ^ 



= ^«(det5)e'^^'=6^ed^ 



-^OceObf + -Oce^b} + -GceObf + ^ce^bf 



K 



K 



K 



+ ^-{deiB)e^'^efUB-')t^,ebf + ^(deti?)6»'"=6,e.(i?-^)f ^e^e + ^K(deti?)e'^^^e^ed^eeee,/ (51) 



Continuing with the expansion. 



13 



d 



;(detS)(<5«<5S - S2St){B-')f^eae + ^^{detB){B-y ^^'^ 



3 
— / 
2 



edfeece^bf 



= --{detB) [{B-')t^ - S'}tr{B-')gj\ e,f - -{detB) [(5"^)^^ - S>{B-')gj 



+^K{detB){B-^^'^ ^ 



Note that the first two terms on the right hand side of (52) are equal to each 
other, and equal modulo a factor of 3/2 to (50). Combining (49) and (52) 
produces a contribution due to the derivative terms of 



2(detS)[4tr(S-)A_(5-)^^_^ 



^bf + -K{detB)iB-^)f-^e''''^edfeeceebf.i53) 



It is convenient to define a new vector field antisymmetric in indices be and 
antisymmetric in indices fe given by 



fe » 

and its single and double traces, given by 



J- = \^^^-^- A = V^ r55) 

^ dx'^f ~ dxbf ' dx ^ dx^f ■ ^ ' 

J ce /e (,/ 

Equations (55) can as well be written in the form 



(56) 

The functional divergence is the divergence in the space of functions at 
each spatial point x in S, which has the interpretation of the next influx or 
outflux of a physical property indigenous to the CDJ matrix elements within 
a 'volume' in the space of functions. The question is what is the source that 
drives this flux relative to the SU{2) directions. One can see that it has to 
do with a certain combination of matter fields and their first-order quantum 
fluctuations encoded in il^. Note that there is a trace only in 5?7(2)_ indices 
in this divergence and not in position label x. Hence the divergence reduces 
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one from nine degrees to freedom to one degrees of freedom per point for 
Eae, and from 81 superficial to one degrees of freedom for eae^hf- 

Thus the total contribution to qi is given, combining the results of (46) 
and (53), by 

2(detB)^efc/ + -K(det5)^(6„ee6/) - (6C^ + 25^trC)ev 

+6(tre)(trC) + 6tr(eC) + GnC^^eaeebf + GQi (57) 

One may in general desire to be able to write this in the form of a func- 
tional covariant derivative, but first let us simplify the linear part of the 
nonderivative terms of (57) involving C in terms of e;,/. 

-(6Cj + 2S'}tTC)ebf + 6(%trC)e6/ + ecjcb/ = 4trCtre (58) 

So, the overall contribution to qi due to quantum gravitational fluctuations 
about the pure Kodama state is given by 

2((detB)^ + 26bftrC)ebf + -/^((det5)^ + 4C»)) (e^e^;,/) + Gn^ 

(59) 

Equation (59) is the total gravitational functional divergence. Note that 
each term acts as a covariant derivative along a particular direction in func- 
tion space, using = AfBl as a connection. Recall that this is the same 
connection utilized in the SU{2) representation of the Gauss' law constraint 
[11]. So the connection gives the rule for parallel transport of the CDJ 
matrix elements both in functional and in position space, two mutually or- 
thogonal sets of directions in the trivializable base space (i;,r) of the fibre 
bundle structure. It will be convenient to use V for the covariant form of 
the functional divergence operator d in order to place these derivatives on 
an equal footing as spatial gradients, highlighting the trivializability of the 
bundle. In this notation we have, for the total contribution to gi, 

(?i = {2{deiB)d''f + 4SltTC)ebf + ^{{detB)dll + 4Cff)eaeebf + GQi. (60) 



5 Expansion of qo relative to the pure Kodama 
state 

This is the only algebraic constraint to be solved, apart from the diffeo- 
morphism constraint, and is nonlinear. We will need to make use of a few 
identities. Starting with the Ansatz 
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= -{^Sae + eae) (61) 



it suffices to calculate the SU{2) invariants of the CDJ matrix. We first 
calculate the trace and variance 



tr^' = -(^ +tre^ 
Far* = (tr^')^ - tr*^ 



- + tre - (-Sab + eab)i-Sba + ^ba) 
+ -tre + (tre)^ -{^ + -tre + tre^ 
= ^ + -tre + (tre)'' - tre"' 

K 

6 4 

— > Var^ = Vare + ^ + -tre (62) 

K"^ K 

The SU{2) trace and variance have a statistical interpretation. 

(i) If the values of the CDJ matrix were correlated to a normal distribution, 

the SU{2) trace would correspond to the mean and the SU{2) variance 
would correspond to the variance of statistical fluctuations about the mean. 
But what is the interpretation of statistical fluctuations in SU{2)? 

(ii) We interpret these to be actual quantum gravitational fluctuations of 
order k cx A, due to the presence of matter fields. These statistical moments 
would not exist for the pure Kodama state for which eab = 0. The pure 
Kodama state can therefore can be viewed in terms of normal statistical 
distributions as a spike of zero width. Therefore we expect that the presence 
of matter in a universe with a cosmological constant should induce quantum 
effects which should be distinguishable by our model. 

(iii) There is a large number of different matter models and combinations 
thereof that can be examined. But whatever the model, the induced quan- 
tum effects should in principle fix the CDJ matrix elements. 

(iv) There is another statistical measure dete, the third-moment, which ap- 
pears in our model. In statistical language, it is a measure of the skewness 
of the distribution. So we need to determine what the analog for skewness is 
for SU{2). Neverthless, it is a quantum gravitational effect which is carried 
over into the scmiclassical level, albeit suppresed by order relative to the 
'statistical mean' and is present due to the cosmological constant. 

This brings us to the final invariant that we will need in the expansion 
of the Hamiltonian constraint about the pure Kodama state, namely the 
determinant 

det^' = ^eabcedef^ad'^be^cf (63) 
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Substituting in the Ansatz, 



-6det^' = €abcedef{-Sad + Ead) {-^be + Efee) (-^c/ + Cc/) 

K/ Ki K/ 



— 1^ "^^abc^defSadhe^cf + eabc^def^adhe^cf + 3k ^ '^abc'^def^ad^be^cf 

+^abc^def^ad^be^cf = 6k"^ + 3K-2(e"^^e„fe/)ec/ + Sk"^ (5^5^ - S}S^)€be€cf 

+eabcedef<^adebeecf = 6k~^ + 6K~^tre + 3K~^Vare + 6dete(64) 
Putting it all together, 



qo = detB(Adet* + Var^) + GQq = detB(6Kdet* + Var^) + GOq 
= detB^K( - 6k~^ - 6K~^tre - Sn'^Vare - 6dete^ 

+detB(6K~^ + ^tre + Vare^ + 
= detB(^-2K~hie - 2Vare - Gndeiej + 0^0. (65) 

The semiclassical part of the Hamiltonian constraint, then, is 

qo = GOo - 2det5(^^ + Vart + 3/tdete) = GQq - detfij^^tre + 2Vare + Adete) (66) 

Note that the mass dimensions in (66) balance, which is a good check on 
the algebraic steps. Rearranging (66) we obtain 

A A^ GA 

tre + -Vare + — dete = j^^o- (67) 



5.1 Expansion of matter contribution relative to "^Kod 

We will assume that the scalar potential F((/>) can be included as a con- 
tribution to the cosmological term A. Since the Klein-Gordon scalar field 
contribution to the quantized Hamiltonian constraint is quadratic in con- 
jugate momenta we should expect a contribution both to ^Iq and to fli. 
Starting with the classical form of this contribution. 

Recalling the results from [11], defining di<pdj(l) = Tij, we have for 
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_ i diT 1 d ,6 . _ i diT 1 d , , 

which has the interpretation of an inhomogeneous quantum matter effect 
proportional to dn/dcf), plus a contribution linear in €ae- 
For CIq, we have 



^^0 = Y + ^Tij^ae^afBiBj 

^1 '6 6 

= Y + ^(BBfTij + eae{\T,^BlBi) + ^T.j BlB}eaeeaf 

6 1 

= Too + -^Tae^ae + -^SabTefeaeebf (70) 

which has the interpretation of an inhomogeneous term corresponding to Tqo, 
the time-time component of the matter energy momentum tensor coupled 
to the background metric {BBy^, plus a part linear in eoe, and a quadratic 
correction. 



6 Full expansion of the constraints 

We are now ready to expand the quantized constraints into a nonlinear 
system of nine equations in nine unknowns. There wil be six equations, all 
linear, from the kinematic constraints and three equations (go = 9i = Q'2 = 0) 
from the Hamiltonian constraint. 

6.1 Total contribution to 

Rewriting (67) maintaining all indices explicit. 



A 

GA^ GA^ ^ 

= ^"^00 + -^TaeGae + -^OabTefeaeCbf (71) 
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The idea is then to transfer all terms linear in Sae to the left hand side of 
(71), maintaining all inhomogeneous and nonlinear terms on the right hand 
side. This yields 

= ^^"^00 — -^[OaeObf — OafObe " -^OabTefjeae^bf — l^^abc^def^ad^be^cf 

(72) 

In (72) ry"^ is not to be confused with the spacetime Minkowski metric r]^'^ , 
since the ae indices take their value in SU{2)-. Rather, it should be taken as 
a kind of S'f/(2)_ metric which corresponds to DeSitter spacetime, corrected 
by the >S'J7(2)_-projected energy momentum tensor of matter. 



6.2 Total contribution to qi 

(II) Moving on to the qi term and dividing through by a factor of 2, 

((detS)5- + 2daetvC)eae + ^((detS)^:^ + 4C$)..ee,; + ^^ + jr.,e^eae = 0(73) 
We now transfer all terms linear in Cae to the left hand side of (73) to obtain 

V'febf = [\B\{6bfiB-')t - iB-')%f + %r,f)^ + 26bftrC]ebf = \^- |v:;e„e6,/(74) 
where in (74) we have made the definitions 

v:i=5:f+4c«) (75) 

and also 

Cff = 51C) - SlC'j. (76) 

We have in (74) transferred all linear terms to the left hand side, leaving the 
inhomogeneous (quantum) term and nonlinear terms on the right hand side. 

The inhomogeneous quantum matter term q serves as the driving source for 
functional divergence relative to the pure Kodama state. 

(III) Lastly, the q2 term reads 

A«^e„e = (77) 

with defined as in (37). Equation (77) is a linear operator on the CD J 
matrix elements. 
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7 Matrix representation of the constraints 



We are now ready to put the constraint equations into matrix form. First, 
let us rewrite the equations (29), (31), (72), (74), (77) for completeness. 
They read, leaving the hats off the operators for simplicity. 



^ead^ae 
G'feae = GQd = 



diT A. 



A"^eae = 0(78) 



the set (78) can be written in the form 



Oltecd = GQab + U- (79) 

where G is Newton's gravitational constant. 

In (79) O is a linear transformation acting on the CDJ deviation matrix, 
which can be represented by a nine by nine matrix of differential operators. 
These differential operators can be thought of as vector fields with respect to 
directions in the base space BoKod of the fibre bundle Eaxod- In a sense, the 
operator O acts as the analogue of a structure group for the fibre bundle, 
which transforms the components of the fibre (section) e^b in a manner 
prescribed by the dynamics of quantum gravity as encoded in the quantized 
constraints. Schematically, 



EcKod = (S, r, e,e(S, r)) (E, r, e„b(S, r), o^^(r; as, 5^)) (80) 
where we have adopted the notation 

0:i{T;d^,dA) ^ Og[r,At; ^, ^) (81) 

for X* = {x,y,z) {t^,t'^,t^) in E and {Af,(f)°') in T. In (81) we have 

highlighted the model-specific matter dependence of the operator O. The 
matrix representations in (79) are given, first for the 'structure' transforma- 
tion matrix O, by 
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It will be important to keep note of the mass dimensions of the components 
of the operator O as displayed above. This will be of significance when we 
invert the operators. The mass dimensions are given by 



[1]=0; [G^1 = 3; h""] = 0; [V""] = 3; [A""] = 0; M = -2 (82) 



yah 



The mass dimensions must balance on both sides of (79). 
The CDJ deviation matrix is given by 

/ ^21 \ 
€32 
ei3 

ei2 

^ab = ^23 
€31 

en 

£22 

V £33 / 

Note how the matrix elements eae are arranged such that the correct linear 
transformation O is implemented. The inhomogeneous term Qab is given by 
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The quantities roo are the S'f7(2)_ analogue of the time-space part of the 
matter energy momentum tensor in dimcnsionlcss units, given by 
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The quantity Qa is the SU{2)- charge, and is of mass dimension [Qa] = 3. 
It is zero for the Klein-Gordon field, but we had put it in for explicative 
purposes. The quantity q is given by 

and has mass dimension [q] =3. This is a quantum effect which manifests 
itself already at the linearized level of the constraints. 

All sources associated with the Einstein's classical metric relativity, namely 
the dynamic components of the stress-energy tensor Tq^, have been made di- 
mensionlcss. However, the sources indigenous to classical Ashtckar variables 
and quantized matter effects, Qa and q are not. This imbalance needs to be 
rectified, since the elements of O^f mix the components of Qab together in 
(79). It will be rectified during the process of inversion of the matrix O. 

Our method of solution to (79) will be by iteration. First we solve 
the linear inhomogeneous equation with the error absent and then we 
compute the error. This error contributes to the inhomogeneous source term 
Qab for the next stage of the iteration, etc. In the process we should build 
a power series in the source Qab and then establish criteria for convergence 
of the series. The inversion of the matrix O appears problematic since its 
elements consist of noncommuting differential operators. So one requirement 
for the existence of solutions to the quantum constraints is that the inverse 
be well-defined. 



7.1 The error vector 

Before we attempt inversion of O^f, let us write down the error vector (,ab- 
The error vector encodes the nonlinearities inherent in the quantized Hamil- 
tonian constraint and forms the driving mechanism by which the CDJ devi- 
ation tensor should converge, under the flow induced by the linear operator 
O, to a unique solution. A lot about the error vector, which is model specific, 
hopefully may be used to assess quantities such as fixed points and stability 
or instability of the solution. The orbit or flow of eae within the fibre of the 
fibre bundle EcKod can be interpreted as a kind of renormalization of the 
SU (2) metric ^ab relative to DeSitter spacetime. The error vector is given 

by 
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In tensor notation this is given by 

^ab = <5al'^6lCl[eae, + 5a2426[eae, (85) 

where ^1 and are given by 

= - [kP^'^eaeebf + A^E-'^-^eade,eecf) ; 6 = —K'f^aeebf (86) 
where we have defined 

/"^'^ = ^ (<5ae5v - Saf5he " J<5a6Te/) ; = le«^^e<^^^ (87) 

The mass dimensions are given by 

[6] = -2; [6] = 1 (88) 

Note that ^2 is independent of the model and is model specific, scaled 
relative to its model-free counterpart by terms of order Tae- 



8 Inversion of O^f 

One strategy to invert the nine by nine matrix O'^ is to split O into diagonal 
and off-diagonal parts O = D + R, with D diagonal and the remainder R 
off-diagonal. We will assume that the diagonal part is invertible since there 
is no commutation of differential operators required. Since the matrix O 
consists of differential operators, either with respect to S or F, its inverse 
will consist of path-ordered integration along directions in the base space 
BoKod = (SjT) of the fibre bundle. Inversion of the diagonal part is tan- 
tamount to finding a basis with respect to which the fibre bundle structure 
(S, r, ^'a(,(E, r)) is manifestly trivialized. That this is possible is a direct 
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consequence of the commutativity of the Gauss' law constraint with the 
quantized Hamiltonian constraint [9], [10], 

[Ga,iJ]=0. (89) 

We then compose the inverted diagonal part with the nondiagonal part 
to form a composite operator [/, which acts on the right hand side. We 
express the result formally as a geometric series of operators U = D~^R, 
which can be put in path-ordered or functional-ordered form in analogy to 
the results of [11] regarding the Gauss' law constraint. Schematically, we 
manipulate (79) into the form, 

Oe = GQ + ^(e). (90) 
First we solve the linear inhomogeneous part of (90) 

Oe^ S — ^ e ~ O-^S. (91) 

where S = GQ represents the source vector consisting of matter-dependent 
terms. We do this schematically via the steps, omitting the hats for sim- 
plicity. 



C = {D + R)e = S — ^De = -Re + S 
e = -{D-^R)e + D'^S — ^ (1 + U)e = D'^S = Q' 

oo 

e = (1 + Ur^Q' = (5^(-l)"C/") Q' = [Pe-^] Q' (92) 

n=l 

where the diagonal part is given by 
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and the off-diagonal remainder R is given by 
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In order to solve the linear part of the constraint equations it is necessary 
only to invert five differential operators, 



= Diag{-1, -1, -1, iG^)'^ (G^f^ {Gsf\v'\^22, A33)"' 
= Diag{-1, -1, -1, (Gf (Gf)-\ (Gf)-\ ^, {^22)-', (Aaa)"^) (93) 

where r]^^ is a c-number. The five differential operators correspond to 
five linearly independent directions within the base space BoKod = (S^T), 
namely three directions in 3-space S (due to Gauss' law) and two direc- 
tions in functional space T (due to qi and q2)- These directions form a 
five-dimensional subspace of BqkocL- Any boundary conditions due to inte- 
gration in these directions should be included as part of the inversion process 
(93), but we will in practice take them to be zero, since we are comparing 
^GKod to "^Kod, the latter for which e^e = V {x, (x)). 

The inversion of D implements two main effects: 
(i) It transforms Qah-, which contains some dimensionful components as- 
sociated with Gauss' law charges and quantum terms of the Hamiltonian 
constraint, into Q'^^ which is dimensionless. This is given by Q' = D~^Q 
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At this stage all elements of the source term of (78), being dimensionless, 
have been placed on equal footing. A convenient physical interpretation is 
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that q and Qa, external structures with respect to Einstein's classical metric 
gravity are initially unsynchronized with respect to Tq^j in BcKod- Nonlocal 
transformations of Qa and q, respectively with respect to S and F, are re- 
quired to bring them into synchronicity. In this sense the energy Too has been 
brought into synchronicity as well, albeit through a local transformation, 
(ii) The inversion of D also determines [/, which can be seen as a normalized 
version of D, the rows normalized by their respective diagonal elements. 
U = U^^ acts on all components of the 'normalized' source Q'^^ identically. 
In effect, Einstein's metric relativity, Ashtekar variables and quantum theory 
have all been placed on the same footing. The matrix U is given by 

-hchd 
—^GcSed 

{Gf)-^Gt-5,25d^ 
{Gfr^Gf - 6c3Sdi 

(A33)-lA-'^ - 5e35d3 / 

The mass dimensions at the linearized level of the flow equation 

e„5 = GC/„t(QLi + ea(,) (94) 

are given by 

[e„e] = -2; K6] = 0; [Q J = 0. (95) 
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8.1 Path-ordered integration in the unified base space 

Before we attempt a general solution of the flow equation let us note a 
general property of the solution. The inversion of differential operators 

with respect to the five orthogonal directions in BQXod corresponds to path- 
ordered integration with respect to these directions, which traces out a path 
in a five-dimensional subspace. 

Let us unify the base space of the bundle into the generalized form 
(S, r) X) where x = Xab is the variable of integration in question, whether 
it be one of the three spatial coordinates or one of the two components of 
the connection A". Then associated with each Xab is a Green's function with 
respect to propagation in position space for each function for all functions, or 
propagation in functional space for a fixed position for each position. Let 
be the set of variables of differentiation comprising the matrix operator i?, 
one variable of differentiation for each nontrivial clement of R, each forming 



26 



a vector field in BQXod- The net effect of R upon tfie source Q' is to mix 
their differentiated components together in a specific way. There are a total 
of 5 X 9 — 5 = 40 such variables i/, but only 3 + 3x3 independent variables 
per spatial point in the full theory (three from space S and nine from A^. 
Therefore the vector fields form an overcomplete basis. 

The general solution to the inversion of D then is given, assuming trivial 
boundary conditions (which would correspond to the pure Kodama state, 
with a fixed point Eae = 0), by 

{D-^USM] = r" dx'abKabiXab, Xab)Scd{Xab), (96) 

where Pab = Pab{^) is a function for which the 'functional' boundary con- 
dition on eae is determined. This is analogous to choosing the boundary, 
in usual differential equations, upon which the boundary conditions on the 
function being solved for are specified in order to produce a unique solu- 
tion. In the field-theoretical analogue the choice can be made judiciously, 
for example such as to impose certain desired conditions on ^cKod- 

The two-index object Kab in (96) is the propagator or 'generalized' 
Green's function or functional heat kernel associated with inversion of the 
component oiD in the direction oixab- Incidentally, the action oiU = D~^R 
on a source vector can be interpreted schematically in the following way, if 
we suppress all dependence upon any spectator variables Oi in Bqi^qji. R cor- 
responds to a generator of transformations of the source term S = S[x, J^] in 
the v directions and D^^ corresponds to the inverse transformation in the x 
directions which are in general different. The effect of the noncommutativity 
of these transformations is to change the u dependence by a x-dependent 
amount, or to mix the coordinates. The path-ordered exponentiation of 
this transformation can be seen as a kind of 'generalized' translation of the 
source within Bxod- Schematically, 

[Pe-^] S[x, ~ S[x, V + r?(x, v\ «)] (97) 

for some nonlocal function t). Hence the effect of (97) is to drag the source 
S from one position in BcKod to another with each iteration, tracing out an 
orbit in BQ^od- The unexponetiated action is given, written in components, 

by 

^ dx'abKabiXab, X'ab)RlWff]Scd{x'ab, ^fe) 

(98) 

We will now write down a general solution for the constraints of general 
relativity. To recapitulate, the constraints can be written relative to a nine- 
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dimensional vector space in some irreducible basis. The fibre of EqkocLi 
namely the CDJ deviation matrix e^e can be written as the components of 
a 9-vector in this space 

£;= (E,r(S),e„;,(r),S)) {x,A'i{x),r{x),eab{x,At,r))- (99) 



9 Iterative general solution for e^e 

The full set of quantized constraints can be written in the form 

0':^e,a = GQ',, + U (100) 

The correction terms ^ab are suppressed relative to the linear parts by 
small numerical factors. Moving on to the construction of the solution, it 
will be convenient first to express the error vector in tensor notation 

^ab = SalSbl^ll + Sa2Sb2^22 + Sa3Sb3^33 (101) 

The straightforward method for solving the constraint equations is to 
solve the linear part first. To linear order one has 

Oabi^cd){0) = G{Qab){0) ^ {^ab){0) = Ct^ab (Qcd)(0) • (102) 

Note in (102) that we have taken initial error (^ab)(o) = 0. We also have 
taken this invertibility as given, since we have from the previous section an 
explicit expansion of U to all orders, as well as its interpretation as a kind 
of nonlinear group transformation. The CDJ matrix to this order is given 

by 

^ab = -l{Sab + l{GA)U:tQ'cd) (103) 

Note that all terms in brackets in (103) are dimcnsionless. The coupling 
constant GA is also dimensionless in units for which h = 1, though it does 
not contain a factor of h. 

Next we compute the error tensor, using the solution for e from (102), 
along with the matter input Q to yield {^ab){o) = ^ab[GU~^Q]. Note that 
this expression also takes into consideration the explicit dependence of Q 
upon Eae- This is given, using (85), (86), and (87) by 
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fK'^ae '^bf J^'^ae '^bf ^ af J^a'e'^b' f' 

where in (104) we have made the identifications 



fTaibia2b2 XX raebf fTa\b\fja2h-2 

'-'gh ~ "gl'^hl^ <^ae '^bf 

-5o26H2{v%{Kl''Utf') + K^'Ubt^tf), (105) 

which is an operator that acts in the appropriate fashion on the Q' vectors, 
and hkewise 

fja,bra2b2a3b3 ^ ^^_^^^^^abcdef (ja^b^fja2b2-(jap ^^Qg^ 

Also, we have omitted the subscripts on (Qa&)(o) i^i (104) and (105) for sim- 
phcity. Taking the error vector (104) for this stage, which is now exphcitly 
a function of the matter source variables, write the new constraint as 



Oa6(ecd)(l) = G(Q'„;,)(o) + (Ca6)(l) [GUQ'^q-^] 

{eab\i) = GUg{Q'^\o) = GKtQcd + Kbi^cdh) [GUQy (107) 

whereupon a solution is obtained for the second stage of iteration. The 
corresponding CDJ matrix is given by 

= -| (<5a6 + liGmgQ'^ + l{GAfu!^u;},''''''''Q'^,bMa2b2 

+liGA)'u£u;t''''''''''QUQUQU) (108) 

The process is then repeated via the identification 

(ea6)(n+l) = G't^ab (<3cd)(n+l) = {(^ab){n+l) + Uah{£.cd){n) [{^ab){n)] (109) 

and the resulting series analyzed either for convergence as well as for stabil- 
ity and nonperturbative behavior relative to the fixed point ^^od foi' which 
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{^ab)n = Vra. At each stage in the iteration procedure (109) one is solving 
a linear equation which can be expressed explicilty in terms of {Q'^i,)o- The 
interpretation is that all the quantities making up the dimensionless vector 
Q' are renormalized. This includes the time-space part of the energy mo- 
mentum tensor of matter Tq^ and the SU{2)_ charge Qa as well as q. The 
space-space part Tij is not renormalized since it is not part of the source, 
however it does affect the renormalization of the source. 



9.1 Sophisticated interpretation 

Now we move on to the sophisticated interpretation. The idea is to define 

a new 9- vector (ab forming the left hand side of (100). The interpretation 
is that if Q'^ represents the unrenormalized source, then (^ab is the fully 
renormalized counterpart, given by 

Cab = 0^tecd = GQ',,+U- (110) 

One then expresses the right hand side of (110) in terms of (^ab- A similar 
iterative procedure is then performed, except this time with a convenient 
physical interpretation. First we must calculate the error vector. We will 
illustrate this to first-order, for simplicity, in order to illustrate the basic 
idea. We will maintain all indices explicit in this simplified example. 
The procedure is as follows: 

(i) We must first eliminate in favor of Cab, which requires the inversion 

eab = U^U- (111) 

(ii) We must now express the components of the error vector explicitly in 
terms of C,ab- Let us redefine them, from (86) 

A2 

Cl = —^Ibf^ae^bf — l^^abc^def^ad^be^cf (112) 

due to the nonlinearity of the semiclassical part of the Hamiltonian con- 
straint, and 

^2 = ~Vffeaee,f (113) 

We can now compute the necessary terms, (iii) Rewrite the quadratic con- 
tribution to ^1 in the following form 

Ibf^aef^bf = IbWae ^bf Ca'e'Cb'f = Ca'e'Cb'/'- (114) 
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This has the interpretation of a hneax transformation of the deformed SU(2)- 
invariant single contracted epsilon tensor. Of course, it is not invariant under 
the transformation U, since U is more general than SU{2), and also due to 
the matter field contribution. We have made the identification 

Kt = I {S^S'} - d-fSl - ^dabTef) Ui:'ulf (115) 
(iv) Rewrite the cubic contribution to in the following form 



dete = -eabcfdef^adebeecf 
= ^abc^defUad fjlj (a' d' Cb' e' Cc' f = '""^ ^ Ca'd'Cb'e'Cc'/' (US) 

Again, (116) can be viewed as a linear transformation of the double ep- 
silon tensors out of SU{2)^. Taking into account the covariant derivative 
operator, we can write (in full-blown form, for completeness) 



Cab = GQ',, + 5al5n (GA|S|-lj^o[C] - (l/6)Ay,«/CaeC6/ " l^^l'K«<iC6eCc/) 

-5a2h2{{l/2)G^,[C] + |(Vfy + r^^)C„eCv|ll7) 

where we have defined 

The sophisticated form of the constraint has a convenient interpretation. 
Note that (117) is exact to all orders and no approximation has yet been 
made. Think of Qah as a vector in an abstract nine-dimensional space that 
corresponds to the full solution of quantum gravity. The right hand side 
can be thought of as a transformation of this vector to a new point in this 
space. It is considerably more convenient to perform an iterative solution in 
this way, since it is easy to keep track of the n^^ iteration by computer for 
arbitrarily high n. One simply makes note of the orbit of this vector under 
the generalized linear transformation and tests for various properties and 
measures of convergence and stability. 

Let us assume that C,ab has been found to the order desired. Then one 
can construct the CDJ matrix from this. It is given, simply by 

"^ab = -jSab + KtCcd (119) 

Since the CDJ matrix can be thought of as an inverse cosmological constant 
(field-valued 'tensor'), the process of iteration of (117) can be viewed as a 
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renormalization of this 'generalized' cosmological constant: a renormaliza- 
tion due to quantum gravitational effects, which manifest themselves at the 
classical level. 

The general solution can then be written as a power series in the cosmo- 
logical constant, in the form 

„ cxD n 
Cab = J ^{GArU^^'^^^'--^-'- n Q'a.h (120) 
n=l k=l 

C,ab in (120) has the interpretation of a renormalization of the source in an 
infinite series in powers of a dimensionless coupling constant GK. The multi- 
indexed operator U^^ consists of products of tj which translate the respective 
sources acted upon within the generalized manifold comprising the base 
space. Each such effect is suppressed by a factor of GK. The operator also 
contains powers of the spatial part of the matter energy momentum tensor 
Tjj , which is separate from the components contained in the source Too and 
Toj. By the use of the Einsteins equations 

oc G/,, (121) 

the power series (120) can be expressed as a power series of positive powers 
of the Ricci curvature tensor R^n and the curvature scalar R = g^^R^v and 
negative powers of the Ashtekar curvature det-B. The full CD J matrix can 
then be written in the form 

^ab = +5:(GA)-f7-^-^^— n Qa.^.y - j 

n=l fe=l 

(122) 

The compact notation on the right hand side of (122) signifies a new kind of 
generating function or group action for the asymptotic expansion, in which 
the 'generator' u acts on the 'parameter' Q' in a kind of tensor representa- 
tion. The specific form of the n*^ iteration of this generator can be explicitly 
carried out based on the entities specifying the matter model in queston. 
One can then make the identification 

1 " 

In this way u generates all of the propagators and vertices of the expansion. 



10 Wavefunction of the universe: Discussion 
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10.1 Mixed partials condition as a boundary condition on 
the wavefunction of the universe 

The mixed partials condition is a consistency condition on the quantization 
procedure which must be consistent with the solution to the quantum con- 
straints and can be viewed as arising from the commutator of the conjugate 
momenta, as in 



[5l(x),7r«(x)] =OVx. (124) 

Equation (124) is simply the condition that the gravitational and the mat- 
ter momenta, independent dynamical variables, have trivial commutation 
relations with one another. This implies, in the Schrodinger representation, 
that 



Plugging in the CDJ Ansatz into (125) results in the condition 



The semiclassical part of (126) does not lead to anything new, however the 
vanishing of the singular term implies the condition that 



*GKod = (126) 



dAf d(t)^ G "dcj)'^ ^ ' 

whereupon the semiclassical matter momentum tTo is directly related to the 

CDJ deviation matrix One can now integrate each a component of 

(127) in the functional space of fields with respect to to obtain more 
conveniently in CDJ variables^ 




In (128) /«[(/>] plays the role of a constant of functional integration with 
respect to A^ = A^{x) and depends entirely upon the matter fields. There 
are no a-priori restrictions on these functions, but it is convenient to pre- 
scribe it as a sort of boundary condition on quantum gravity with a physical 

^This makes use of the nonconvcntial calculus in [11] in which spatial and functional 
integration and differentiation are commuting processes. 
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interpretation. The function fa can be thought of as the conjugate momen- 
tum for the matter fields in the absence of gravity. One could thus choose fa 
to correspond to the conjugate momentum for the matter field quantized in 
the absence of gravity (limc^+o)- If one takes the governing quantum theory 
to be quantum field theory in flat Minkowski spacetimc, then one can impose 
the observable semiclassical limit for this theory as a boundary condition on 
the gravity-coupled quantum theory.^ The solution to the constraints then 
extrapolates this boundary condition nonlinearly as a of backreaction on 
DeSitter spacetime. 

There are two commuting processes occuring in (128). There is partial 
differentiation with respect to the functional dependence of €ae in the direc- 
tion of the matter variables cf)", as well as integration in the direction of the 
gravitational variables Af. 



10.2 Input from the semiclassical limit into "^cKod 

The canonically determined wavefunction for the generalized Kodama state 
can be written, in the reduced phase space quantization scheme, as the 
exponential of the starting Lagrangian evaluated on the solution to the con- 
straints [11]. 



^GKodiX"', r] = exp dt [{HG)-' d^x^aeBiAf d3^7ra0"](.129) 

When the quantum constraints are solved, the semiclassical-quantum cor- 
respondence is identically satisfied and the wavefunction corresponding to 
reduced phase space and Dirac quantization must be the same^. It is clear, 
in the solution of the constraints for a matter-coupled theory, that both 
the CD J matrix and the semiclassical matter conjugate momentum tTq, 
must contain dependence both upon the gravitational and the matter con- 
figuration space variables. This is a natural consequence of the interaction 
between gravity and matter. 

In order to acquire some intuition on the functional boundary conditions 
on the wavefunction of the universe it is instructive to deduce the form of 
the wavefunction (129) in the limit when there is no interaction. There are 
two extremes: First for pure gravity with cosmological term and no matter 
fields present the wavefunction should reduce to the pure Kodama state. 
This can be most directly seen in CDJ variables from the solution of the 
corresponding Schrodinger equation 



^By the same token, the self-duahty condition can be see as a boundary condition on 
the coupled theory in the absence of matter fields. 

''This is a conclusion drawn by the current author, based on the results of [11] and [12]. 
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A abc ^ ^ ^ 

Q^'^'^^'^' SA^x) 6A%x) SAlix) 



abc 



SAf{x) 6A%x) 



B'cix) 



GKod 



which leads to a nontrivial solution for the chosen operator ordering of mo- 
menta to the left of the coordinates 



^GKoM, = ^Kod[A] = e-se^^)-^^^^. (131) 

The pure Kodama state has semiclassical orbits that correspond to DeSitter 
spacetime [2], [3]. 

The other extreme is the case when there is pure matter in the absence 
of gravity. Then one might expect the wavcfunction of the universe to 
reduce to the ground state of the quantized matter theory in the Schrodinger 
representation in Minkowski spacetime. In this case, the function iTa [A, (/)] 
fa{4>) serves as the conjugate momentum corresponding to the WKB state. 
Take for simplicity the massive Klein-Gordon field in Minkowski spacetime. 
The time evolution of the state is governed in the Schrodinger picture, by 
the functional Schrodinger equation 

ih—^KG[i>] = H[(f>, -iS/S(t>]^KG[4'] = E^KG[<t>]- (132) 
where the Hamiltonian operator is given by 

where E is the energy eigenvalue of the state. Let us set the self-interaction 
potential V^cp) to zero for simplicity. Hamiltonians quadratic in the dy- 
namical variables allow for a simple expression for the quantum mechanical 
wavcfunction '^[4>] for the state in terms of a mode number basis of har- 
monic oscillators. As is known, the ground state in this basis corresponds 
to a Gaussian. The corresponding ground state for '^[(f>] is given by 

^kg[<P] = exp[-^ J d^xd^y(Pix)G{x,y)<P{y)\ (134) 
where the kernel G{x, y) satisfies the relations 

/j3„ pip-{x-y) 
A I (135) 

The energy eigenvalue E is just the trace of the kernel G{x,y), given by 
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E = ll dh: J ^Vl^f + m2, (136) 

the zero-point energy of the scalar field (^(x). The function / in this case is 
determined by acting on the state 



tt{x)-^kg = -i^TTT^'^KG = / d^yG{x,y)4>{y) = /[(/'(x)]*/^g(137) 

5(j){x) IJy, -I 

with a more complicated functional form for / when the potential F(^) 
is nonharmonic. The form can be deduced from the effective action for 
a quantized theory in flat Aiinkowski spacetime, either by some kind of 
expansion of the solution to (133), or via the path integration approach 

"UmatteAm] = e^^/'^)^^// ['^(^)l = J D^e-^^t>\ (138) 

For a harmonic potential V = (l/2)m^^^ the effective action can be com- 
puted expactly from (138) to be a Gaussian, which can then be expressed in 
terms of the boundary values of the classical solution. In the more general 
case of a nonharmonic potential one resorts to a pcrturbative loop expansion 
about a reference configuration of the field 4>{x,t), usually chosen to corre- 
spond to a stationary point of the action S'[(^] = S[(l)ci] = Sci[ct>{T)], where 
T is the time corresponding to the final conditions of the classical solution 
(pel as defined on the final hypersurface 

Z{J) = e^^(-^) = (Det(v2 - m^))-'^''e-^^^^t'('''^^Wp[J d^x V[S/SJ{x)] 

^-jj d^xd^yJ{x)l^{x - 2/) J(y)](139) 



exp 



One then computes the value of / via the identification Te// = —W{J) + 
J J ■ (f), and consequently 



whereupon one makes the identification f{(t>) = 5T(,f f/8(f). This is the semi- 
classical limit forming the input into the quantities Q'ae- This result can 
be applied to a wide class of models more general than the Klein-Gordon 
scalar field. 

When the Klein-Gordon field is coupled to gravity, the corresponding 
Schrodinger equation takes a more complicated form given which can be 
written as 
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The imprints of the scalar field model are rearranged when quantized on 
the same footing as gravity, but one should hope that the solution to (141) 
should reduce to that for (130) in the limit of pure gravity and to that for 
(132) in the limit of pure matter. This provides a physical basis for fixing 
the constant of integration /„ in the mixed partials condition corresponding 
to (141), which is given by [11] 

n{x) = fim) + ^[j^5X-^{x))^^eae[X'f{x),ct>{x)] Vx. (142) 

To make the boundary conditions more apparent one can expand the 
conjugate momentum variables in (129) about the corresponding boundary 
conditions. This is given by 

*ae = - {\5ae + e„e) ; TT^ = fa + ^{J ^^Mf ) ^Cae- (143) 

As long as these expansions (143) mutually consistent, then the resulting 
wavefunction should automatically incorporate the appropriate semiclassical 
limit into the coupled quantum theory. The source terms in the quantum 
constraints arising from have the interpretation of a backreaction on 
the spacetime determined by the pure Kodama state, which is DeSitter 
spacetime. Substituting (143) into (129), one obtains 

Jto Js 

Equation (144) consists of a pure gravitational ground state ^cKod^ a matter 
part which drives statistical quantum fluctuations from this ground state, 
which is given by the wavefunction for the quantized matter field in the 
absence of gravity and a part encoding the backreaction effects as deter- 
mined by the solution to the constraints, which is labeled by the semiclas- 
sical matter momentum / = fa{(t>)- The function forms an input via the 
inhomogeneous source vector in the constraint equations 

Hi = UD4r- Qa = fa{Tarp<l^; J^o = roo(/); ^^i = -^E||^! = 02(/Xl45) 
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10.3 Generalized Kodama state for the full theory 

The generalized Kodama state, modulo prefactors due to path integration, 
can be expressed as [11] 



^GKod[A, = exp / {{HGr'^aeBiA^ + -irq 



M 



(146) 



C„.b=0 



Incorporating the mixed partials condition into the matter momentum [11], 



7r{x) = f[4>{x)] - 



Gd(p 



(147) 



We have switched from d to 5 notation regarding the matter field for consis- 
tency with the functional one-form Sep. The reader should bear in mind that 
all field-theoretical infinities have already been factored out, so S/S(j){x) is re- 
ally a partial derivative with respect to (f) at the point x and not a functional 
derivative. 



'^GKodf [A, 4)]=e 
exp 



-6(?iGA)-i/cs[^] 



exp 



(hG)-^ J (fx J 8X1 



(148) 



So we see that there is a contribution to the gravitational sector from the 
matter field due to the mixed partials condition. To see more clearly the 
nature of this contribution, note that there are two operations occuring at 
each spatial point x. There is differentiation with respect to ^, holding 
fixed, followed by integration with respect to 0, again at fixed A". This 
looks like a minisuperspace effect, but is in fact still the full theory. The 
differentiation and integration over the functional space of (j) are inverse 
operations to each other, and thus cancel each other out. The next result 
is that the matter momentum makes a contribution to the wavefunction 
exactly matching that from the gravitational sector. Equation (146) leads, 
applying 122 to the results of [11], to 



^GKodf [A, 0] = ^Kod[A]e^v [-^ 



exp 



{l2hGK)-^ / (fx / (5X°*^(GA)"C/; 



yaifoia2 b2---anbn 



ab 

n=l k=l 



n«u[A<^,/]ii49) 



There are a few things of note regarding the generalized Kodama state: (i) 
It is labeled by /, a function purely of the matter fields, as a nonseparable 
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basis [11]. The function / can be thought of as the semiclassical matter 
momentum in the absence of gravity, and can be chosen such as to impose 
the proper semiclassical limit of the fully coupled theory; (ii) The gravita- 
tional sector is expressible as an asymptotic expansion, on the final spatial 
hypersurface Sy, relative to the pure Kodania state in powers of GA; (iii) 
the functional integration over the space of fields T occurs at each fixed po- 
sition X on T,T to produce a functional of the fields which is then integrated 
over position d^x; (iv) The operators I[^bia2b2...anbn ^^^^ functional 
dependence of the matter source terms Q'^j^ ~ {Toi, Qa, q) upon the Ashtekar 
connection A'^; (v) The asymptotic expansion is in terms of products of the 
local matter matter SU{2) ® Dif feo® H charges propagated around Bxod, 
with coefHcients involving Tij. This might be preferable, due to the dimen- 
sionless coupling constant, to an expansion in powers of curvature found in 
many higher-derivative theories of gravity. 

10.4 Network analogy to Feynman diagrams 

There is another interpretation available for the series expansion (149). If 
one thinks of the source Q'^j^ = Qai,[x, A, 0] ~ g' as a matter-dependent prop- 
erty smeared over the base space Bxod of Exod and the operator U as a gen- 
eralized propagator and the error as a self-interaction vertex parametrized 
by the coupling constant (GA)", then one can visualize the following process 
in which the source traces out an orbit which can be projected onto the base 
space Bxod- 

The source q' = q'^^^ propagates a disturbance in functional space from 
Xo to xi where it undergoes a scattering at a 3-point and 4-point vertex. 
The resulting disturbance combines with the source already present at xi 
to form a new source q" . This new source propagates from to ^2 where 
it undergoes another scattering and combines with the source present at ^2 
to form the source g'" and the process repeats itself. The number highest of 
source particles at a particular stage proliferates as 3" and the total matter 
source, q"" ' behaves as though it is undergoing a kind of random walk 
through BQKod- The walk is not really random, it is predetermined by the 
dynamics of the quantized constraints. At each stage the source accumulates 
more and more particles. One hopes that the total source to infinite order 
should be convergent. 
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